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INTRODUCTION 


Radomes  significantly  reduce  missile  guidance  accuracy  by  producing  boresight  error 
(BSE),  an  error  in  target  direction  as  determined  by  a  seeker  or  radar  system.  Despite 
decades  of  development  and  use,  radomes  are  a  problem  because  newer  missiles  require 
small  angular  rates  of  change  of  BSE  and  because  newer  seekers  have  greater  frequency 
bandwidths.  The  electromagnetic  design  process,  for  boresight  error  and  bandwidth,  is 
largely  empirical  because  radome  shapes  are  non-separable,  making  analysis 
approximate.  The  need  for  improving  radome  performance  motivates  efforts  to  improve 
analytical  methods. 

The  purpose  of  the  research  is  to  develop  improved  methods  for  analyzing  radome 
electromagnetic  effects.  The  main  goals  are  as  follows: 

HF.  Develop  a  geometric  theory  of  diffraction  (GTD)  for  dielectrics  and  in 
particular  hollow  dielectric  shells 

Extend  the  moment  method  for  use  with  hollow  dielectric  shells. 

The  GTD  is  intended  to  be  relatively  simple,  and  the  moment  method  is  intended  to  be 
accurate  although  more  complex  than  the  GTD. 
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SUMMARY 


Section  1  reviews  the  concept  of  transmittance,  the  complex-valued  wave  distribution  in 
a  radome  bounded  volume,  and  give  measured  and  computed  values  of  transmittance  for 
pointed,  axially-sym metric  radomes  and  for  wedge  radomes.  The  transmittance  has 
intensity  variations  that  are  not  described  by  analytical  methods  that  utilize  a  flat  sheet 
transmittance  to  locally  approximate  the  radome.  These  flat  sheet  methods  omit  wave 
mechanisms,  namely  guided  waves  and  vertex  scattering.  To  quantify  these  mechanisms 
the  moment  method  was  applied  to  a  finite  slab  for  grazing  incidence.  We  postulated 
guided  waves  and  waves  scattered  by  a  slab  edge.  The  amplitudes  and  phases  of  these 
waves  were  determined  by  interpreting  maxima  and  minima  computed  by  the  moment 
method  in  terms  of  the  postulated  waves.  Thus  we  have  solved  numerically  a  canonical 
problem  for  a  geometric  diffraction  theory  of  dielectrics. 

Section  2  extends  the  moment  method  to  a  hollow  cone  with  a  scalar  Green's  function 
approach,  and  Section  3  gives  numerical  results.  Section  4  gives  measured  results  for  a 
hollow  cone.  Section  5  compares  the  numerical  and  computed  results.  Agreement  is 
generally  good.  The  shapes  of  measured  and  computed  curves  agree,  but  the  maximum 
phase  is  somewhat  in  error.  Polarization  dependence  is  correct.  The  formulation 
decomposes  the  cone  into  cells  that  are  sectors  of  annuli.  This  decomposition  utilizes  to 
advantage  the  constant  values  of  fields  in  a  cell. 


Section  6  describes  a  tensor  Green's  function  formulation  for  a  hollow  cone.  The  cells 
are  spheres.  This  decomposition  treats  polarization  more  simply  than  does  the  annular. 
Section  7  compares  measured  and  computed  results  for  a  rod  and  a  ring. 

Appendix  A  shows  the  theoretical  equivalence  of  the  two  formulations. 


During  the  coming  year  we  plan  to  test  further  the  two  formulations  for  off-axis 
incidence  on  hollow  dielectri$  wedges  and  on  hollow  dielectric  cones.  Numerical  data 
will  be  generated  for  cqyer^and  tested  against  measured  values.  The  numerical  data  will 
be  interpreted  to  develop  geometric  diffraction  theory  for  hollow  cones,  which 
approximate  radomes. 
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1.  GEOMETRIC  DIFFRACTION  THEORY  FOR  DIELECTRICS 


1.1  APPROACH 

The  central  concept  is  the  wavefront  distribution,  of  phase  and  intensity,  produced  in  a 
radome-bounded  region  by  an  externally  incident  plane  wave.  This  distribution,  called 
the  transmittance  T,  depends  on  incidence  direction,  wavelength,  and  polarization,  and 
the  radome  configuration  and  composition.  The  transmittance  strongly  influences  the 
output  voltages  of  an  antenna  enclosed  by  a  radome.  Physically,  the  phase  variations  of 
T  are  aberrations;  in  fact,  wavefront  tilts  are,  to  a  first  approximation,  proportional  to 
boresight  error.1  Figure  1-1  suggests  aberrations.  Analytically,  the  transmittance  T  is  a 
function  that  is  a  factor  in  the  integrand  of  a  diffraction  integral  for  the  patterns  of  an 
enclosed  antenna.  The  quantities  that  are  measured  with  an  antenna  are  intensity  |ER|2 
and  phase,  arg  ER,  where  ER  is  a  rectangular  field  component 

ER  =  TE*F  dA, 

where  E1  is  the  incident  field  component,  F  is  the  receiving  near-field  distribution,  and 
the  integration  is  over  the  region  such  that  F  is  non-negligible. 

The  transmittance  T  helps  in  radome  development.  For  externally  incident  waves, 
nearfield  values  of  T  directly  link  boresight  error  (via  wavefront  tilts)  and  the  radome 
shape  and  its  thickness  variations.  Thickness  variations  and  shape  modifications  are 
limited  in  extent,  but  they  are  the  main  parameters  available  for  design.  The  link 
provided  by  T  is  more  direct  than  that  between  boresight  error  and  farfield  patterns 
because  diffraction  smooths  phase  variations.  Of  course,  farfield  patterns  are  signifi¬ 
cant  for  missile  performance;  in  fact,  boresight  error  and  aberrations  were  related 
analytically  through  polynomial  expansions  of  T  in  diffraction  integrals.1 

We  approach  a  GTD  for  T  by  measuring  and  computing  transmittance  near  hollow 
dielectric  shells,  hollow  wedges,  and  a  flat  slab.  A  slab  is  an  idealized  model  of  half  of 
an  axially-sym metric  radome  for  high  incidence  angles.  We  seek  to  decompose  the 
complicated  transmittance  into  a  sum  of  simpler,  constituent  waves,  as  in  Figure  1-2. 
Identifying  constituent  waves  and  estimating  their  magnitudes  and  phases  generates 
conical  solutions,  an  important  part  of  GTD.  The  measurements  are  done  with  a  small, 
scanning  probe  antenna.  The  calculations  are  done  by  two  radome  analysis  methods,  one 
a  direct  ray  method  and  the  other  a  surface  integration;  in  addition,  we  apply  the 
moment  method  fo*  infinitely  long  cylinders.2 

The  approach  to  the  second  goal,  extending  the  moment  method,  is  to  analyze  a  hollow 
cone.  Two  formulations  are  used.  One  uses  a  scalar  Green's  function;  the  other  a  tensor 
Green's  function.  The  two  are  equivalent,  but  the  numerical  procedures  differ.  The 
scalar  Green's  function  method  decomposes  the  cone  into  circular  cylinders  and 
decomposes  the  cylinders  into  angular  sectors.  The  tensor  Green's  function  method 
decomposes  the  cone  into  spheres.  The  formulations  differ  in  two  practical  ways.  One  is 
in  the  number  of  cells  because  of  cell  size.  The  other  difference  involves  polarization 
dependence. 
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1.2  BACKGROUND:  ANALYTICAL  METHODS  BASED  ON  FLAT  SHEET  APPROXI- 
MATION 


A  simple  way  to  analyze  radomes  is  sketched  in  Figure  1-3.  This  direct  ray  method 
associates  a  single  incident  wave  normal  (ray)  with  a  point  on  the  receiving  aperture. 
The  method  is  adequate  for  designing  radomes  big  enough  to  enclose  an  antenna  that  is 
an  order  of  magnitude  larger  than  the  wavelength,  provided  the  aperture  also  is 
considerably  larger  than  the  wavelength.  However  this  method  fails  for  smaller  antennas 
and  radomes;  in  these  cases,  a  surface  integration  approach  has  been  successful.^ 
Figure  1-4  suggests  surface  integration,  which  can  be  interpreted  in  terms  of  rays 
converging  to  an  aperture  point. 

In  both  the  direct  ray  and  surface  integration  methods  the  transmittance  for  each  ray  is 
approximated  by  the  values  for  a  plane  sheet  that  locally  approximates  the  curved 
radome,  as  in  Figure  1-5.  Surface  integration  more  densely  samples  the  variation  of  the 
surface  normal  direction.  Figure  1-6  shows  the  variat.  n  in  the  normal.  For  the 
simplest  case  of  axial  incidence,  a  converging  set  of  rays  spanning  the  tip  has 
transmittance  that  is  a  sum  of  terms  with  distinct  individual  transmittance.  One  term 
arises  from  each  converging  ray  in  Figure  1-4.  The  value  of  each  term  depends  on  the 
direction  of  the  surface  normal  for  the  associated  ray  because  the  normal  direction 
determines  the  angle  between  it  and  the  field  and  thus  the  resolution  of  the  field  into 
parallel  and  perpendicular  components.  The  direct  ray  method  sparsely  samples  the 
variation  of  the  surface  normal  direction,  but  dense  sampling  is  significant  near  the 
radome  tip  where  the  normal  varies  rapidly. 

Note  that  transmittance  depends  on  position  within  the  radome.  Figure  1-6  suggests  how 
this  variation  occurs.  On  the  Y-axis  a  direct  ray  projected  from  the  radome  gives  the 
transmittance  value  for  parallel  polarization,  and  on  the  X-axis  it  gives  the  value  for 
perpendicular  polarization.  The  variation  with  position  requires  surface  integration  to 
describe  accurately  the  mixture  of  terms. 

In  addition  to  position  dependent  polarization  variation  relative  to  the  surface  normal  for 
fixed  polarization  direction,  another  variation  occurs  when  the  field  orientation  is 
changed  from  horizontal  to  vertical.  Consider  Figure  1-7  which  compares  measured 
phase  with  values  computed  by  the  direct  and  surface  integration  methods  for  off-axis- 
incidence.  it  shows  that  surface  integration  is  more  accurate  than  the  direct  ray 
method.  It  also  suggests  how  a  polarization  change  influences  transmittance,  even  for  a 
fixed  point  in  the  volume  bounded  by  the  radome. 

1.3  MORE  BACKGROUND:  EXAMPLES  OF  TRANSMITTANCE  FOR  AXIALLY 
SYMMETRIC  AND  WEDGE  RADOMES 

Although  Figure  1-7  and  the  results  in  Reference  3  suggest  that  surface  integration  is 
more  accurate  than  the  direct  ray  method,  these  results  are  for  radomes  that  were,  to  a 
good  approximation,  a  half  wave  thick.  For  wide  frequency  bands,  surface  integration, 
as  described  in  Figure  1-4,  also  has  limited  accuracy.  To  illustrate  this  point, 
Figures  1-8  and  1-9  show  boresight  error  of  a  radome  with  thickness  0.125  in.  for 
frequencies  from  8  GHz  to  18  GHz.  The  antenna  was  an  array  of  two  horns.  For 
8-12  GHz  each  horn  aperture  was  1.95  in.  by  1.7  in.  The  results  show  rather  small 
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ANMSOA 


Figure  1-3.  The  Direct  Ray  Method.  A  Single  Incident  Wave  Normal  (or  Ray)  is 
Associated  with  a  Point  on  the  Receiving  Aperture.  The  Radome  is 
a  Surface  of  Revolution,  but  only  a  Plane  is  Sketched. 


MHO* 

Figure  1-4.  Surface  Integration  Method.  The  Radome  is  a  Surface  of  Revolution. 
Integration  is  over  a  Portion  of  the  Incident  Wavefront. 
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Figure  1-5.  Approximate  Description  of  Radome  at  a  Point  by  a  Flat  Sheet. 
Multiple  Internal  Reflections  are  Included. 

discrepancies  for  some  frequencies  but  larger  discrepancies  at  other  frequencies.  For 
example  in  Figure  1-8,  for  the  E-plane,  frequency  12  GHz,  and  gimbal  angle  -16°  the 
discrepancy  is  0.3°.  The  discrepancies  for  10  and  14  GHz  are  smaller  for  the  same  angle, 
but  they  are  larger  for  8  GHz. 

To  understand  the  discrepancies  in  boresight  error,  we  have  compared  measured  and 
computed  transmittance  values  for  the  radome  of  Figures  1-8  and  1-9.  Figure  1-10 
shows  an  example,  for  12  GHz,  16°  gimbal  angle,  and  the  E-plane.  The  point  of  this 
figure  is  that  it  illustrates  the  magnitudes  of  transmittance  discrepancies  (in 
Figure  1-10)  that  correspond  to  those  of  boresight  error  in  Figure  1-8. 

We  also  considered  another  radome,  with  thickness  tapered  between  0.45  in.  and  0.40  in. 
This  thickness  is  a  half-wave  wall  at  10  GHz,  but  it  is  20  percent  greater  at  12  GHz. 
Calculated  and  measured  transmittance  values  for  an  axial  position  that  would  be 
occupied  by  an  aperture  are  shown  in  Figure  1-11.  Consider  two  horns  centered  at 
X  =  ±A.  The  phase  difference  between  these  two  points  is  greater  in  Figure  1-11  than  in 
Figure  1-10.  Therefore  boresight  error  discrepancies  should  be  larger  than  for  the  case 
of  Figure  1-10,  and  boresight  error  measurements  have  verified  this  expectation. 

To  study  this  radome  (of  Figure  1-11)  further,  we  measured  the  transmittance  through 
the  radome  enclosed  volume.  Figure  1-12  shows  power  transmittance  (intensity)  for  16° 
gimbal  angle. 

This  graph  shows  transmittance  variations  between  +1  dB  and  -4  dB.  An  antenna  would 
be  illuminated  by  a  wave  that  is  non-uniform. 
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ure  1-6.  The  Variation  in  Direction  of  the  Surface  Normal  N  Produces  Transmittance  That 
Depends  on  Position.  Even  for  a  Fixed  Incident  Polarization,  the  Transmittance  Varies 
with  the  Normal  Direction  from  Parallel  to  Perpendicular. 
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Figure  1-9. 


i 


9 


Figure  1-10.  Phase  Delay  and  Power  Transmittance  for  Pointed  Radome,  Length  to 
Diameter  Ratio  1.8,  Thickness:  0.125  in.  Frequency:  12  GHz;  E-plane;  Gimbal 

Angle:  16°.  Measured  ( - );  Direct  Ray  Calculation  (-  -);  Surface 

Integration  with  49  Points  Spaced  by  0.75"  (x).  Surface 
Integration  with  49  Points  Spaced  by  1.0"  (°). 
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Figure  1-11.  Phase  Delay  and  Power  Transmittance  for  Pointed  Radome,  Length  to 
Diameter  Ratio  1.8.  Thickness  Tapered  from  0.40  in.  to  0.45  in.  Gimbal 
Angle:  16°,  Frequency:  12  GHz.  Measured  ( — ).  Surface 
Integration  with  49  Points  Spaced  by  0.75  in.  (x);  with 
1  in.  spacing  (•). 
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Figure  1-12.  Power  Transmittance  at  12  GHz  for  a  Radome  with  Dielectric 
Constant  2.6.  Thickness  0.40"  to  0.45",  for  the  Field  Polarized  in 
the  X-Y  Plane.  Each  Curve  is  Labelled  for  a  Particular  Axial 
Position  of  the  Probe  Aperature  and  Unit  Transmittance  Value 
(0  dB)  is  at  the  Corresponding  Axial  Position. 
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In  Figure  1-12  the  transmittance  near  the  wall,  for  Y  <  0,  is  shown  in  Figure  1-13.  This 
figure  shows  a  quantity  As,  where 

As  =  2  7r  (kg-k0sin)-1.  (1-1) 

This  quantity  is  the  spacing  of  fringes  produced  when  a  plane,  free-space  wave  interferes 
with  a  coherent  slab-guided  wave.  Therefore  the  transmittance  in  Figure  1-12  appears 
to  result  to  some  extent  from  a  guided  wave,  as  was  suggested  in  Figure  1-2.  The 
quantitative  question  of  guided  wave  magnitude  remains  as  does  the  reason  for  the  -4  dB 
transmittance  minima.  We  return  to  these  questions  subsequently. 

We  see  that  the  transmittance  of  a  radome  is  complicated.  The  direct  ray  and  surface 
integration  methods  omit  some  propagation  mechanisms  suggested  in  Figure  1-2.  We 
have  seen  in  Figure  1-13  evidence  for  guided  waves.  Reference  4  gave  evidence  for 
guided  waves  in  two  and  three  dimensional  radomes.  In  an  earlier  report,  Reference  5, 
we  determined  guided  wave  amplitude  on  a  finite  slab  by  postulating  functional  forms  for 
guided  waves. 

To  study  the  effects  of  polarization  and  dielectric  constant,  we  measured  the  intensity 
transmittance  of  a  ceramic  radome,  dielectric  constant  5.5,  for  the  H-plane;  see 
Figure  1-14.  The  minima  Me  very  deep  so  an  antenna  would  be  immersed  in  a  very  non- 
uniform  field.  However  the  minima  spacing  are  irregular  and  do  not  help  identify  a 
guided  wave. 


Figure  1-13.  Transmittance,  from  Figure  1-12  at  End  Points  of  Scans  % 
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To  separate  the  effects  of  curvature  we  measured  the  intensity  behind  a  hollow  wedge. 
Figure  1-15  shows  |T|2  for  a  curved  wall  wedge  for  the  E-plane.  Because  T  is  high  for 
parallel  polarization,  the  minima  are  shallower  than  those  in  Figure  1-14. 

Moment  method  calculations  were  made  for  a  wedge  with  flat  walls;  see  Figure  1-16. 
The  deep  minima  and  their  spacing  suggest  the  exponential  decay  (with  distance  from  the 
wall)  that  characterizes  slab-guided  waves.  In  Figure  1-16  the  decay  with  distance  is 
approximately  that  for  a  slab,  but  only  near  the  wall. 

We  conclude  that  guided  waves  exist  on  radomes.  Their  significance  depends  on  the 
thickness-to-wavelength  ratio  because  this  ratio  largely  determines  the  transmittance  of 
the  direct  wave.  For  a  half-wave  thickness,  the  direct  wave  transmittance  magnitude  is 
large;  but  at  a  frequency  20  percent  higher  than  that  for  a  half  wave  the  direct  wave 
transmittance  is  low.  The  perturbation  caused  by  the  guided  wave  depends  on  its  own 
magnitude  as  well  as  that  of  the  direct  wave. 

However,  the  guided  and  direct  waves  do  not  seem  to  be  an  adequate  description. 
Guided  waves  are  exponentially  dumped  so  they  fall  to  small  magnitudes  within 
approximately  a  half-wavelength  from  the  wall.  Therefore,  in  Figure  1-12,  guided  waves 
do  not  seem  to  explain  the  transmittance  minima  and  maxima  at  distances  approximately 
a  wavelength  from  the  radome.  An  additional  mechanism  seems  necessary. 

A  plausible  mechanism  is  a  vertex-scattered  wave.  Some  experimental  evidence  exists 
for  a  vertex  wave  in  wedge  radomes.®  However,  the  situation  is  less  clear  for  an  axially- 
symmetric,  pointed  radome.  We  reduced  transmittance  discrepancies  by  adding  a  vertex 
wave  to  direct  ray  calculations7,  but  this  result  was  obtained  before  surface  integration 
was  developed.  The  vertex  wave  may  be  a  construct  that  corrects  the  direct  ray 
method. 

Because  the  transmittance  minima  and  maxima  in  Figure  1-12  seem  significant  we 
explore  a  vertex  wave  as  a  possible  explanation.  The  approach  is  to  analyze  a  flat  finite 
slab  with  the  moment  method. 

1.4  A  GEOMETRIC  THEORY  OF  DIFFRACTION  FOR  FINITE  DIELECTRIC  SLABS 

This  section  considers  grazing  incidence  on  a  finite  dielectric  slab  in  order  to  model 
approximately  some  of  the  waves  scattered  by  a  hollow  shell.  The  model  is  partial 
because  the  field  in  a  shell  contains  a  constituent  wave  that  propagates  through  the 
illuminated  side,  but  this  wave  does  not  occur  for  a  slab  at  grazing  incidence.  Some  or 
all  of  the  other  waves  shown  in  Figure  1-2  may  be  present.  In  particular  a  guided  wave  is 
excited  in  either  the  slab  or  the  shell;  of  course,  the  magnitude  of  the  guided  waves 
differ.  The  slab  is  an  idealization  also  because  it  has  not  circumferential  curvature; 
consequently,  neither  the  incidence  angle  nor  the  polarization  direction  vary.  Moreover, 
the  discontinuity  at  the  slab  edge  differs  from  that  at  a  radome  tip,  but  the  discontinuity 
distinguishes  both  from  an  infinite  slab.  The  slab  approximates  one  side  of  an  axially 
symmetric  radome,  and  a  wedge  approximates  both  sides. 

The  purpose  of  this  section  is  to  explain  power  transmittance  minima  like  those  in 
Figure  1-12.  Because  of  these  minima  a  seeker  antenna  would  be  illuminated  by  a  non- 
uniform  field  so  boresight  error  and  antenna  patterns  would  be  affected.  The  minima  are 
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not  given  accurately  by  surface  integration  calculations;  see  Figure  1-11.  Several 
mechanisms  may  exist.  Reflection  from  the  left  side  of  the  profile  in  Figure  1-12  is 
ruled  out  by  ray  diagrams.  Although  Figure  1-13  shows  evidence  for  guided  waves,  the 
lowest  order  slab  guided  mode  would  decay  to  negligible  values  within  a  wavelength  of 
the  radome  wall.*  Therefore,  another  mechanism  seems  necessary  to  explain  the 
minima,  which  are  separated  from  the  wall  by  more  than  a  wavelength. 

To  study  possible  mechanisms  we  start  by  computing  the  field  in  and  near  a  two- 
dimensional  slab  with  the  moment  method  of  Richmond.  Coordinates  are  shown  in 
Figure  1-17.  The  slab  has  infinite  extent  in  the  plane  orthogonal  to  the  figure,  and 
incidence  is  in  the  X-Y  plane.  The  slab  is  approximated  by  a  set  of  infinite,  circular 
cylinders,  and  the  fields  at  the  centers  of  the  cylinders  are  determined.  Figure  1-18 
shows  the  total  field  intensity  at  each  cylinder  center  and  the  total  field  outside  the 
slab. 

To  check  the  computed  field,  measurements  were  made  with  a  dipole  probe.  The  dipole 
scanned  several  paths  parallel  to  the  slab.  Measured  intensity  is  shown  in  Figure  1-19. 
Agreement  is  good.  In  Figure  1-18,  the  fringe  spacing  for  Y  =  0.54  in.  agrees  with  As 
computed  from  Equation  1.  This  spacing  is  evidence  that  guided  waves  were  excited  on 
the  slab.  The  shallow  oscillations  with  period  0.5  in.  result  from  interference  of  guided 
waves  travelling  in  opposite  directions. 

*A  calculation  later  in  this  section  gives  the  decrease  with  distance  from  the  wall. 


Figure  1-17.  Coordinates  for  a  Slab 
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Figure  1-19.  Measured  Intensity  Near  Slab  (Length  10  in.;  Thickness  0.25 
Wavelength  1.26  in.;  Dielectric  Constant  2.6). 


The  scattered  field  magnitude  (the  total  field  minus  the  incident)  is  shown  in 
Figure  1-20.  The  field  is  not  uniform,  indicating  that  the  scattered  field  in  the  slab 
contains  more  constituents  than  just  a  guided  wave.  Since  the  spacing  of  peaks  in 
Figure  1-20  is  As  to  a  good  approximation,  the  field  in  the  slab  apparently  consists 
mainly  of  a  guided  wave  and  a  plane  wave  with  free  space  propagation  constant.  In 
addition,  a  small  amplitude  reflected  guided  wave  is  also  present.  The  free  space 
propagation  constant  k0  rather  than  nkQ,  when  n  is  the  index  of  refraction,  is 
remarkable.  Apparently  Icq  occurs  because  the  slab  is  thin,  not  a  volume  in  which  n 
modifies  the  wavelength. 

Tc  test  the  guided  wave  hypothesis  and  to  estimate  its  excitation,  we  examine  the  curve 
fo.*  Y  =  0  in  Figure  1-20.  Consider  the  oscillation  of  |ES|.  The  value  varies  from  2.7  to 
approximately  0.1,  where  |E*|  has  unit  value.  Now  smooth  the  variation  with  period 
0.5  in.;  |ES|  varies  from  2.3  to  0.1,  or  1.2  ±1.1.  We  assume  the  amplitude  1.2  represents 
that  of  a  guided  wave  propagating  in  the  positive  X  direction.  Now  consider  the 
variation  of  a  slab-guided  wave  with  Y.  In  the  slab,  the  field  has  the  form 

E„.  =  1.2  cos  u  y;  (1-2) 

Gi 

outside  for  Y  0  the  field  has  the  form 

E  =  0.96e“v(y'b)  (1-3) 

Go 

where  v  is  computed  from  boundary  conditions  as  1.12  inches  and  b  is  half  the  slai. 
thickness.  The  scalar  amplitude  0.96  was  computed  from  Equation  1-2,  with 
u  =  1.65  inches.  The  values  computed  from  Equations  1-2  and  1-3  are  shown  in 
Figure  1-21.  For  comparison,  estimates  of  guided  wave  amplitude  derived  from 
Figure  1-20  and  additional  calculations  by  the  moment  method  are  shown  in  Figure  1-21. 
The  agreement  in  Figure  1-21  is  additional  evidence  for  a  slab  guided  wave  and  verifies 
the  estimate  of  excitation  magnitude. 

The  apparent  value  of  k0  rather  than  nk0  discourages  a  simple  ray  picture  that  assumes 
nkQ  for  an  internal  ray;  however,  outside  the  slab  a  free-space  plane  wave  (plus  an 
evanescent  guided  wave)  seem  reasonable. 

To  determine  the  form  of  the  wave  that  combines  with  the  guided  wave,  we  graph  the 
phase  increment  of  the  scattered  field  between  adjacent  points  on  the  X-axis.  Inside  the 
slab  at  Y  =  0,  Figure  1-22  shows  wide  variations.  Outside,  the  free  phase  increment  is 
reached  approximately  for  Y  =  0.94  in.  The  phase  increment  k0Ax  is  71.5°  for 
AX  =  0.25  in.  At  even  larger  distances  (6  in.)  the  amplitude  of  the  scattered  field  is 
below  0.1. 

To  obtain  a  comprehensive  view  of  the  fields  near  the  slab  we  plotted  the  contours  of 
computed  phase  and  intensity  of  the  scattered  field  for  grazing  incidence.  See 
Figures  1-23  and  1-24.  The  total  field  phase  and  intensity  contours  are  plotted  in 
Figures  1-25  and  1-26. 


1-21 


INCIDENT 


NORMAL 


AML051  X  (CM) 

Figure  1-23.  Computed  Scattered  Field  Phase  Contours  Near  1/4  in.  Thick, 

10  in.  Wide  Slab  with  k=  2.6;  Wavelength:  1.26".  Incidence:  Grazing. 

The  Contours  are  Spaced  by  360°. 

The  phase  of  the  scattered  field  shown  in  Figure  1-23  suggests  a  cylindrical  wave  from 
the  edge  nearest  the  source.  Within  the  slab,  the  phase  changes  rapidly  near  X  equal 
5  in.  This  behavior  was  shown  graphically  in  Figure  1-22;  see  the  graph  for  Y  =  0  near 
X  =  5  in. 

We  conclude  that  the  scattered  field  consists  of  guided  waves,  a  cylindrical  wave  from 
an  edge,  and  an  approximately  plane  wave  at  larger  distances  from  the  slab.  However 
the  analysis  does  not  explain  why  the  second  constituent  of  Es  should  travel  as  xQ  x  in 
the  slab. 

Consider  the  total  field  in  Figure  1-26.  For  the  region  outside  the  slab,  we  assume  the 
total  field  consists  of  the  incident  field  and  a  cylindrical  wave  centered  at  the  origin, 
which  is  at  X  =  0,  Y  =  0.  The  field  is  thus 
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where  r  was  approximated  by  x  +  (y2/2x),  and  C  is  as  yet  an  undetermined  constant.  To 
test  Equation  1-4  we  computed  loci  such  that  the  two  terms  are  in  phase.  Figure  1-27 
shows  two  such  loci  superimposed  on  moment  method  calculations.  To  obtain  this 
agreement  we  choose  C  to  be  -1.  That  is,  the  cylindrical  wave  suffers  a  phase  reversal 
on  its  excitation.  Since  amplitudes  were  not  computed,  the  magnitude  of  C  was  not 
determined  by  using  Equation  1-4. 
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Fig  <re  1-24.  Computed  Intensity  in  Field  Scattered  by  Slab.  The  Intensity 
Levels  are  Compressed  in  the  Lower  Picture. 
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Figure  1-25.  Phase  of  Total  Field  for  Slab 


In  Figure  1-26,  the  maxima  and  the  minimum  near  the  slab  resemble  the  variations  in 
Figure  1-12.  Therefore  the  analysis  of  the  slab  seems  relevant  to  analysis  of  radomes  so 
that  some  additional  analysis  on  estimating  amplitudes  near  the  slab  seems  warranted. 

Now  consider  Figure  1-28.  It  shows  the  variation  of  the  fields  with  Y  the  transverse  co¬ 
ordinate  for  X  fixed  at  5  in.  For  the  total  field  E  computed  by  the  moment  method, 
maxima  occur  at  Y  equal  2.4  in.  and  4.8  in.,  and  a  minimum  is  Y  equal  0.6  in.  The  two 
maxima  in  |E|  coincide  with  bright  regions  in  Figure  1-26,  and  the  arrows  indicate 
maxima  positions  determined  by  assuming  a  cyclindrical  wave,  see  Equation  1-4. 

For  the  scattered  field,  |ES|  computed  by  the  moment  method  is  also  shown  in 
Figure  1-28.  We  see  a  minimum  for  Y  equal  0,  despite  the  maximum  of  |E|  of  a  guided 
wave  and  a  wave  with  propagation  constant  *0. 

The  magnitude  of  the  scattered  field's  guided  wave  constituent  |Eq  |  was  graphed  in 
Figure  1-21.  It  is  shown  in  Figure  1-28  on  another  scale. 

The  magnitude  of  the  cylindrical  wave  constituent  |Eq|  is  shown  also  in  Figure  1-28. 
The  value  1.1  for  Y  equal  0  was  estimated  from  the  interference  behavior  of  the 
scattered  field;  Figure  1-20.  For  values  of  Y  greater  than  0.6,  |Eq  |was  taken  as  |ES| 
because  Eq  has  decayed  to  negligible  values.  For  Y  equal  0.6  we  assumed 
|EC|  =  |es|  -  | Eg |  because  the  value  of  |ES|  was  nearly  zero. 

The  phase  variation  at  Y  equal  5  in.  is  shown  in  Figure  1-29.  The  moment  method  value 
of  arg  Es  is  well  matched  by  the  assumed  cylindrical  wave  phase  except  at  Y  =  0,  where 
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Figure  1-27.  Computed  Loci  of  Constructive  Interference 

the  guided  wave  amplitude  is  high.  This  guided  wave  component  is  apparently  present  in 
the  moment  method  results,  but  it  is  of  course  absent  from  the  cylindrical  wave 
constituent  of  the  scattered  field. 

The  preceding  analysis  leads  to  a  geometric  diffraction  theory  description  of  scattered 
fields  near  a  finite  width  slab  for  grazing  incidence.  This  description  is  summarized  in 
Figure  1-30.  The  field  consists  of  the  waves  in  Table  1-1. 

Table  1-1.  Constituent  Waves 


WAVE  TYPE 

VARIATION 

MAGNITUDE  AT  X 

EQUAL  5  IN. 

Incident 

exp  (i*0x) 

1 

Guided  wave  (forward) 

exp  (i  Kgx ) 

1.2 

Guided  wave  (reflected) 

exp  (-i«gx) 

0.2 

Cylindrical  wave 

exp  (*0r  +  *) 

1.1  for  Y  =  0 

0.37  for  Y  =  2  in. 

We  thus  have  solved  a  canonical  problem. 
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Figure  1-28.  Fields  Near  a  Slab  for  Y  Equal  5  in.  E  is  the  Total  Field,  and 
Es  is  the  Scattered  Field,  Both  Computed  by  the  Moment  Method. 

Eq  is  a  Guided  Wave  Constituent  of  the  Scattered  Field  and  E^ 
is  the  Cylindrical  Wave  Constituent.  Frequency  9.4  GHz. 
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Figure  1-29.  Computed  Phases  of  Scattered  Fields.  Moment  Method  (•); 
for  Cylindrical  Wave  with  ir  Added  (+) 
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Figure  1-30.  Constituent  Fields  Near  a  Slab  for  Grazing  Incidence. 
The  Guided  Waves  have  Propagation  Constant  kg.  The  Cylindrical 
Wave  Has  Propagation  Constant  of  Free  Space  kj,. 

These  Waves  are  Geometric  Diffraction 
Corrections  to  Direct  Waves. 


2.  SCALAR  GREEN'S  FUNCTION  FORMULATION 


2.1  INTEGRATION  OF  MAXWELL'S  EQUATIONS 


The  starting  point  is  Maxwell's  equations.8  The  symbols  have  their  usual  meanings,  and 
the  time  variation  has  the  form  exp  (— icot);  thus, 

VxZ  =  iwuH  (2_1) 

o~ 

VxH  =  -iucE  (2‘2) 

A  scattering  object  is  represented  by  considering  the  permeability  €  to  be  a  function  of 
position.  To  solve  Eqs.  2-1  and  2-2,  let  the  total  fields  E  and  H  be  decomposed  into 
incident  fields,  without  the  scatterer,  and  the  scattered  fields;  thus, 

E  =  E1  +  E5  (2-3) 

H  =  2I  +  HS  (2_4) 

where  the  I  superscript  denotes  incident,  and  the  s  denotes  scattered. 

Note  that 

Vxji1  =  iw^H1  (2"5) 

VxH1  =  -iue  E1  (2"6) 

—  o~ 

Now  rewrite  Equation  2-2,  with  e  =  e  +  eQ  -  eQ,  subtract  Equation  2-5  from  Eq.  2-1  and 
Eq.  2-6  from  Eq.  2-2,  and  obtain 

„  .•  ws  (2-7) 

Vx£  =  iu>nQti 

VxHS  =  -iwji  ES  -  iw(€  -  c  )E  ^2  ^ 

—  o  o 

Eqs.  2-7  and  2-8  show  that  Es  and  Hs  satisfy  the  free  space  form  of  Maxwell's  equations, 
except  for  the  term  -ia>(  €  -  cq)E,  which  is  a  polarization  current 


J  =  -ioi€  ^c-l)E 
o  ^ 


2-1 


where  <  is  the  dielectric  constant.  If  we  apply  the  equation  of  continuity,  the 
polarization  charge  is 


P=  -  «  v.  <M)£ 
o 


(2-10) 


5 

A  wave  equation  for  g  results  by  applying  the  curl  operation  to  Equation  2-7  and 
eliminating  H  with  Equation  2-8.  Thus 


Vx  V  x  ES  -  k 
~  o 


(x-l)E  =  ia/MJ 


(2-11) 


Notice  that  the  left  side  contains  £s  but  the  right  side  has  £.  Since  Eq.  2-11  can  be 
expressed  with  J  on  the  right  side  to  replace  we  expect  to  derive  ES  from  potentials  as 


ES  =  -  +  ioj^. 

where  the  potentials  are 

$  =  (4?rO  lJp  g  dv 
A  =  (pQ/4ir)  J"  J  gdv 

and  g  is  the  free  space,  scalar  Green's  function  r  1  expikr,  with  r  the  distance  between 
integration  and  observation  points. 

Although  the  field  in  Eq.  2-12  seems  intuitive  with  the  currents  in  Eqs.  2-13  and  2-14, 
the  classical  deviations9  start  with  wave  equations  (like  Equation  2-11)  that  contain  a 
single  field  quantity  such  as  Es,  not  Es  and  E*.  The  use  of  J  suppresses  this  distinction, 
which  may  be  questionable  when  differential  operators  are  used  to  derive  Equations  2-13 
and  2-14  from  Equation  2-11.  Therefore  we  give  a  direct  integration  of  Equation  2-11. 
The  procedure  follows  that  of  Stratton  and  Chu10. 

s  *  1 

In  Eq.  2-11  let  £  represent  E  ,  Q  =  where  &  is  g  or  r  expikr,  and  a  is  an  arbitrary 
constant  vector.  Apply  Green's  identity. 


(2-12) 

(2-13) 

(2-14) 


J"  (Q-VxVx  P  -  P-VxVx§)  dV  =  J"  (PxVxQ  -  §xVxP)*n  dA 

V  s 

The  left  side  can  be  written  as 


(2-15) 


a  •  [kj  J  (ic-l)<eEdV  -  V'  J  <p V  •  <M)  EdvJ  -  j  (n  •  ES)  (a  •  V<P)  dA  =  L.S. 


(2-16) 


where  the  integration  is  bounded  by  a  large  sphere  andgby  a  small  sphere  about  the 
observation  point.  The  derivation  used  the  facts  that  V  •  E  is  V  •  (*-l)E  and  the  identity 

?fa .  w)  =  V  •  [(a  •  W)ES  j  -  (a  .  v<p)  V  ES 

For  the  right  side  we  obtain 


R.S.  =  JL  *  J*  ^  x  (Es  xn)  +  i<j)n0  (n  x  HS)J  dA 
S 

From  Equations  2-15,  2-16  and  2-17,  we  obtain 


(2-17) 


k2  J  <k-i)  <p<SV  -  V'  J*  <P V  ■  (K-i)  EdV  =  J  [  (n  •  ES)  V<p  +  V<px  (ESm)  {2_18) 
V  V  s 

+  iu>n  v  (nxHS)  dA 
o  "  J 

g 

The  right  side  of  Eq.  2-18  tends  to  4  tt  £  when  integrations  are  done  over  a  large  sphere 
and  over  a  small  sphere  about  the  observation  point.  The  integral  over  the  large  sphere 
tends  to  zero  because  of  the  properties  of  g  and  the  fields,  and  the  integration  over  the 
small  sphere  gives  4  -n  £  .  We  obtain 


4  TT  Es  =  k2  JVi)  g  £dV-  v'Jgv  •  [(K-l)  e]  dV  (2"19) 

which  is  the  result  obtained  by  assuming  Equations  2-9,  and  2-12  through  2-14  as  a 
starting  point. 

The  integral  equation  for  the  electric  field  is  obtained  from  the  definition  of  the 
scattered  field  in  Equation  2-3. 

The  integral  equation  is 


=  EI 

where  Es  is  given  by  Equation  2-19. 
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SIMULTANEOUS  ALGEBRAIC  EQUATIONS 


The  integral  equation  is  solved  for  the  interior  of  the  dielectric  object  by  changing  it  to 
a  set  of  simultaneous  algebraic  equations.  A  hollow  cylinder  is  divided  into  rings  and 
subdivided  into  cells  as  in  Figures  2-1,  2-2  and  2-3.  A  hollow  cone  is  subdivided  into 
rings  as  in  Figures  2-4  and  2-5,  and  the  rings  are  divided  as  in  Figure  2-3.  The  integral 
equation  is  evaluated  at  the  center  of  each  cell  to  generate  as  many  equations  as  there 
are  cells.  At  the  center  of  a  particular  cell,  labelled  with  index  m,  the  integral  equation 
for  each  rectangular  component  is 

_  s  _  I 

mm  m  (2-20) 


for  lSnuSN,  where  N  is  the  number  of  cells.  4  is  a  sum  of  contributions  from  all  cells, 


ES  =  E1 
mn  m 


This  equation  is  put  into  matrix  form  by  writing 


E  -  S  (zS  /E  )e  =  E1  , 
m  Z.  \  mn  nj  n  m' 

Sc  E  =  E1 
mn  n  m 


where 


c  =  1  -  (es  /e  ) , 

mm  '  mm  m/ 


(2-21) 


(2-22) 


(2-23) 


(2-24) 


C  =  -  /E 
mn  mn  n  (2-25) 

The  next  step  is  to  evaluate  the  scattered  fields  E^,n*  Formulas  are  derived  in  the 
following  sections.  Diagonal  and  off-diagonal  matrix  elements  are  separately  given,  and 
the  vector  and  scalar  potentials  are  separately  evaluated. 
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Figure  2-1.  Dielectric  Ring  and  Coordinate  System 
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Figure  2-2.  Subdivision  of  Hollow  Cylinder  Into  Rings.  An  observation  point  has 
coordinates  (x*,  y't  z');  it  may  be  inside  or  outside  the  cylinder. 
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?.3  SCATTERED  FIELD  FOR  DIAGONAL  ELEMENTS 

2.3.1  VECTOR  POTENTIAL 

The  scattered  field  from  the  vector  potential  is 

ES  =  k2  (k-1)  (4ff)  7  g  E  dV  .  (2-26) 

Let  us  assume  horizontal  polarization  both  incident  and  within  the  ring.  This  assumption 
is  reasonable  for  thin  rings  and  axial  incidence.  For  the  cell  labelled  m  we  have 


=  ^  <*_1)  <4,r)’1  f  *  Ev  ^ 

xm  J  x  m 


To  simplify  notation  let 


I.  =  fgE  dV 
Am  J  x  m 

Since  the  integration  extends  only  over  the  cell  labelled  m. 


f  2  2  2*7” 

Lp  "  pm  ~  2ppm  cos  '  V  V  J 


(2-27) 


(2-28) 


(2-29) 


with  integration  limits  Pm  -6 p  <  P  <  Pm  +  8p  ,  zm  -  Sz  <  z  <  z™  +  8z, 

<  0  <  An  +  ^  Introduce  abbreviations  6  =  -  An>  4  =  z  -  zm,  and 

a  =  p  -  pm.  Then  on  assuming $-  is  small,  we  have 


r 


r 

m 


+  pp 


m 


(2-30) 


Am 


6z  6<p  <5p  ikr 

j  J  j  Ex  rm  e  m  (p^  +  cr)  cos  B  d<7  dd  d£ 

-5z  -<5<^>  -<5  p 


(2-31) 


so  that 


-  - ^ 


The  (actor  cos  0  in  the  integrand  arises  from  first  expressing  the  x-component  of  field  in 
polar  coordinates,  expanding  cos  (0  +  4>m)  and  sin  (0+  4>m),  and  then  using  small  angle 
approximations  for  sin  0  and  cos  0.  Now  assume  rm  is  small;  that  is,  <r  and  £  are  small. 
Moreover,  assume  Ex  is  constant  Em  in  the  cell.  Thus  we  can  expand  the  exponential,  to 
obtain 


T  -  T(1>  T(2>  • 

I .  -  I,  -  I .  ; 

Am  Am  Am 


(2-32) 


where 


(2) 

I\  ’  =  ikE 
Am 


m 


6z  6  <6 

J  S' 

-<5z  -6<p  -6p 


6p 

/  » 


m 


a)  cos  9  do  dd  df 


(2-33) 


This  integral  is,  if  we  omit  cr  in  pm  +  cr , 

(2) 

=  ik  (8  5z  op  p  6^)  E  (2-34) 

Am  m  m 

The  factor  in  parentheses  is  the  volume  of  the  cell. 

The  term  1^  is  more  complicated.  We  omit  limits,  which  equal  those  in  Equation  2-31 
Am  --1/2 
r. ;  = 

Am 


=  ///  [(pm  *  c)  %  * 


(0  ~  cr)  cos  9  dcr  dfl  d£  E  (2-35) 

m  m 


2  2  2*1 
1  *  cr  *  t  J 

The  integral  over  £  leads  to  a  logarithmic  function  so 

rAm  =  //  (P  m  ~  a>  C0S  9  l0g  dd  E 

(sWz)]1/2,  and 

(pm  +  cr)  Pm  02J1^2.  Assuming  s+  «5z2,  to  expand  the  square  root, 


(2-36) 


where  R  is 


where  s+  is 
we  obtain 


1  + 
o-2  ♦ 


lAm  *  ff  ^ pm  +  008  9  *°g  'L  (s  /25z)2]}  ^  d&  Em  (2-37) 

*  f  f  (p  +  cr)  cos  9  log  (26/s  )2  dcr  d0  E 
J  J  m  +  m 

+  f  f  (s  /26z)“  (p  +  cr)  cos  9  do  d0  E 
J  J  *  m  m 


We  assume  that  cos  0  l  and  cr«Pm  to  carry  out  the  integrals.  We  finally  have  that 
the  x  component  of  the  scattered  field  is 


4m  ‘ 


“W  2 


(2-38) 


2-9 


where 


2  2  2  2 
s2  =  6p  +  pm 

Equations  (2-27),  (2-32),  (2-34)  and  (2-38)  give  the  scattered  field  in  the  cell  labelled  m 


eL  -  “2  <«-1> 


i  —  k  6z  6p  p  5<b  +• 
t  'm  ' 


pm 


12tt6z 


o  5p<54>  s‘  --J  [l  -  log  (s9/26z)"J 


.(2-39) 


m 


For  brevity  we  omit  the  details  of  the  integrations  over  <r  and  0.  The  integral  was  done 
by  integrating  over  local  polar  coordinates.  In  effect  then  our  cells  are  small  right 
circular  -cylinders  not  the  truncated  wedges  defined  by  pmd0d<r.  This  approximation 
procedure  is  described  near  the  end  of  Paragraph  2.3. 


2.3.2  SCALAR  POTENTIAL 

The  scattered  field  from  the  scalar  potential  is 


where 


~s 


i-a 


(2-40) 


I2  =  fg  V  •  (k-1)  E  dV  .  (2-41) 

To  represent  the  dielectric  constant  we  use  the  unit  step  functions  as  functions  of  p 
only.  That  is,  we  shall  understand, 


(k-1)  =  (K-1)  [H  (Pm  -  op)  -  H  (Pm  -  5p)] 


(2-42) 


where  on  the  right  side  (k-1)  is  a  constant,  but  on  the  left  k  is  a  function  of  p.  Because 
we  have  assumed  Ez  =  0,  it  is  unnecessary  to  specify  the  z  dependence.  Finally,  k  is 
assumed  independent  of  <t>.  Therefore 
ikr 

I2  »  J  ■£—  V  •  (k-1)  fH  (Pm  -  dp)  -  H  (Pm  +  <5p)J  E  (1^  cos  <p  -  1^  sin  <p)  dV  (2-43) 
where  1  and  ^  are  unit  vectors  in  the  r  and  directions  respectively. 

We  assume  div  E  =  0,  in  the  cell  so 


I 


2 


(*-i)  E 

Cl 


(2-44) 
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where 


*  /  g6  (p  -  tp)  cos  4>  dV  ,  (2-45) 

-  J  m  m 

li2)  »  /g<5  (p  -  6p)  cos  £  dV  ,  (2-46) 

£  y  m  m 

=  f  g  cos  <p  dp  d0  dz  .  (2-47) 

The  first  two  integrals,  1^  and  $),  approximately  canceL  To  see  this  consider  iM. 
Assume  kr  is  small  and  expand  expiicr.  Thus 

-  f  (r  1  +  ik)  6  (p  -  <5p)  cos  <J>  dV 

i.  j  m  m 


The  second  term  of  the  integrand,  with  factor  ik,  gives  -4ik  (pm  -  5p)  5z  S4>  cos  4>m. 
This  term  is  omitted  because  it  approximately  cancels  the  term  -4ik  (pm  +  Sp)  Sz  8<t> 
that  arises  from  The  first  term  with  factor  r*1,  gives 


r(D 


-  -1/2 

cos  <t>_  1 1 «/p„)-  +  («p/p„r  *  (i 


m 


m' 


(2-48) 

cos  6  dO  df  ; 


recall  9  and  £  were  defined  earlier.  Following  integration  on  C  we  have 


r(U 


where 


and 


fl--^jcos^  [(fin  -  In  P  )  cos  9  d9  dt 
'  ^zn '  J 


*Sz  lom-JP>9 


m 


2q 


\p 


2  A  2 

<5p 


The  evaluation  of  1^  proceeds  in  the  same  way.  The  result  is  that 


42)  ■ 


where 


^l+^£^cos^tt  JlfLn  Q+  -  In  Q_)  cos  9  d9  df 


Q, 


(a  +■  fip)  6“ 

(q  *  &z)  ^  {pm  p} 


m 


w 


(2-49) 


(2-50) 


(2-51) 


(2-52) 


(2-53) 
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Now  consider 


where 


f’-f’  ■  (x-^) 


dd 


(2-54) 


2.2  p  2  -  5p2 

L  -  3  ~5z  ,  g2  ■  m _  a4.lz5P. 2 

±  2  2 


m 


4s 


e  ± - r  e  . 

S  P-rn 
m 


(2-55) 


Inspection  of  L+/L-  shows  that  it  tends  to  unit  value  for  small  9.  For  example,  at  9  -  0 
the  value  is  1;  for  9  =  0.2,  the  value  is  1.046.  Therefore  the  logarithm  is  small  so  the 
contributions  of  ilP  and  1$)  from  E®  are  omitted. 


Am 


To  evaluate  $*),  expand  the  exponential,  expikr,  in  the  integrand  as  1  +  ikr.  The  first 
term  is  omitted  because  it  is  independent  of  the  observation  point  coordinates  so  the  V' 
operator  yields  zero.  From  the  second  term  we  obtain 


-  Jr1  cos  <p  dp  d<p  dz 


(2-56) 


Again  we  use  substitution  z  -  zm  =  (,  P  -  Pm  =  =  6  to  obtain 


~l/2 


1^  =  cos  <b  f  f  f  cr^  +  (p  +  <t)  p  9~“  +  cos  9  dp  d<?  d£  .  (2-57) 

2  m  J  J  J  m  m  -* 


We  assume  9  is  small  so  cos  9  « 1.  Integration  on  C  gives 


43)  »  cos  (Pm  II  log 


1  + 


2  2 

«v  +  8>  pm  9  ** 

46  I2 


4  6 


r  2  2* 

J  (p  +  9)  p  9  +  a 
pm  m 


do-  d# (2-58) 


Because  wall  thickness  is  much  smaller  than  radius,  we  omit  the  term 


<j  p  9  /4<5  z 
in 


in  the  integrands.  The  integrations  on  cr  and  5  give 


T(3) 


7T  p  COS  <b 

m  m 


.  2  2  <5p  6<£  (6Z)““  2  (2-59) 

log  (s  /26  z)  s  -  — — -  s  cos  <p 

m  J  m  3  m  m 


where 


2 

s 

m 


K  ♦'„>*). 


To  obtain  the  field  apply  the  V'  operator  to  1^).  The  p  and  components  are  as  follows: 


+ 


(2-60) 


+ 


P 


dp  6<b 

m 

2 

37r6z 


sin  <b  E 
m  xm 


(2-61) 


From  Equations  2-60  and  2-61,  the  contribution  of  the  scalar  potential  to  the  x 
component  of  the  scattered  field  is 


z  <*>  =  -  fi-v 

xm 


m 


2  A\ 

cos  © 


m 


(2-62) 


(-gr)  w  + 


xm 


where 


and 


(2-63) 


,  .  2 
f  =  p  op  6<p/{ 3ir  oz) 

m 


(2-64) 


2.4  SCATTERED  FIELD  FOR  OFF  DIAGONAL  ELEMENTS 

2.4.1  VECTOR  POTENTIAL 
The  scattered  field  is 

Es  =  k2  (k-1)  (4  it)-1  J  gEdV  (2-65) 

where  the  integral  extends  over  the  shell,  and  g  is  r-1  expikr,  where  r  is  the  distance 
between  integration  and  observation  points.  The  shell  is  subdivided  into  cells,  as  shown 
in  Figures  2-1  and  2-3;  the  subdivision  is  into  rings,  defined  by  increments  along  the 
z-axis  of  symmetry,  and  further  division  of  each  ring  into  angular  sectors.  The  center  of 
a  cell  has  coordinates  (pm,  0m,  zm)  for  an  internal  observation  point;  integration  cells 
have  coordinates  (Pn,  0n,  zn).  In  these  coordinates 

r2  =  P2  +  P2m  -  2  ppm  cos(0  -  0m)  +  (z  -  zm)2.  (2-66) 

We  let  p  =  pn  +  cr,  z  =  zn  +  £,  and  0  =  0n  +  6,  and  utilize  a  binomial  expansion  so  that 

r  =  B  +  o-  (  ]  B"1  +  £(zn  -  zm)B-1,  (2-67) 

where  B 2  =  pj  +  p2m  -  2pnm  Cnm  +  (zn  -  zm)2- 

and  [  ]  =  Pn-  Pm  cos  (0n  -  0n). 

For  horizontal  polarization  and  axial  incidence  the  volume  integral 

I  =  fg  ExdV  (2-68) 

becomes 

I  =  2n  eikB  B-l  PnEnfeikcr  l']  /Beik^zn  *  Zm)/B  d  <rdi  (28<t>)  (2-69) 

where  the  integral  over  <f>  yielded  280.  We  also  approximated  the  factor  Pn  +tr which 
multiplies  the  exponentials  by  pn.  This  approximation  means  that  radome  thickness  is 
smaller  than  radius.  We  obtain 

E^  =  k2  (k-1)  (2/7 r)  5p808z  2  En  PneikBB-l  sine  k  (  ]  8pB“l  sine  jkZnmSzB"1  (2_70> 

where  znm  is  zn  -  zm,  and  sinez  is  z_1  sinz. 

2.4.2  SCALAR  POTENTIAL 

The  scalar  potential  generates  a  scattered  field 

E®  =  (4tt)-1  V’  Io  (2-71) 

s  4 


2-14 


where 


I2  =  / g  V.(K-l)£dV,  (2-72) 

and  V1  is  taken  at  the  observation  point.  The  integral  extends  over  the  cells  described  in 
the  preceding  section,  so 


h  =  %J  gV-(«-l)IndV.  (2-73) 

The  function  (k-1)  is  constant  within  the  dielectric  and  thus  in  every  cell;  it  depends  only 
on  radius,  which  changes  from  one  cell  to  another.  We  describe  the  inner  and  outer  radii 
by  step  functions  so  that 


(k-1)  =  |k-1|  [H(pn-  5p)-H(Pn  +  Sp)], 


(2-74) 


where  H(x)  is  zero  for  x  negative  and  unity  for  x  positive,  and  pn  is  the  radius  of  the  cell 
with  index  n.  For  horizontal  polarization  we  obtain 


h  =  f  g  V-  |k~1|  [H(pn-3p)-H(pn+  6p)j  Eon(c^l^-s^l^|  dV, 


(2-75) 


where  1^  and  lp  are  unit  vectors.  The  symbols  |k-1|  denote  the  magnitude  of  the 
dielectric  constant  within  the  dielectric. 


Within  the  dielectric  cells  we  assume  V*E  is  zero  because  the  field  is  assumed  constant. 
The  derivative  on  p  of  the  divergence  operator  leads  to  Dirac  6  functions,  but  the 
derivative  on  <P  operates  only  on  sin  <t>, 

because  k  -1  is  independent  of  <f>.  We  obtain 

I2=  Ml  £EnI2n,  (2-76) 


where  the  integral  over  any  cell  is 


(2-77) 


where 


iW 

2n 

I<2) 

2n 

I<3) 

2n 


JgS(pn  -  <5p)  cosApd  p  dd>dz; 
fgS(Pn  +  dp)  cos  4>pd  pdpdz; 
J'g  cos d  pd<t> dz. 


(2-78) 

(2-79) 

(2-80) 
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As  before 


r  =  [p 2  +  p2m  -  2pp  mcos(0  -  <*>m)  +  (z  -  zm)2] 
with  z  -  zn  =  £  ,  and  P  -  p  n  =  o- ;  a  binomial  expansion  gives 
r  =  s  +  bs-1cr  +  znms_1£ 

where  b  is  Pn  -  Pm  cos«>n  -  and  z^  is  zn  -  zm. 

To  evaluate  lg)  let  <t>  be  <f>n  +  9  ,  where  e  is  small  in  each  cell,  so  cos(0n  - 
approximately  cosc^.  The  integrals  over  6  and  are  well  known,  so  we  have 

1^  =  85p  S<^6z  e^^^cos^n  sine  (kb  Sps-1)  sine  (kznm8zs-1)» 
or  approximately 

1^  =  88p8<f>  Sz  eikss_1  cos^ 

For  iM,  we  have 

l(l)  _  e*kss  1  J"  eikb(o-/s)eikznm(  £/s)  8(p-+  gp)  cos(<f>n  +  0)d0  df 
=  4  8  zS4>  (pn  -8p)  cos4>n  e*kss-1  e-^  8P^S, 

where  sine  (k  znm8z/s)  has  been  approximated  by  unit  value. 

In  addition, 

*2ri  =  4  8  z  80  (pn  +  8p)  cos  4>n  e‘kss_1  e*kb  Sp/s. 

Therefore 

j(l)  _  j(2)  -  -8  §  z  8<f> cos  0ne‘kss-1  [Spcos  (kb  8p/s)  +  i  p  nsin  (kb8ps-1]  » 

2n  2n 

with  small  argument  approximations,  we  obtain 

1^  -  1^  =  -88z  84>8p  cosc^ne*kss"^  (1  +  ikpn  bs-1) 

The  electric  field  from  I^j),  Equation  2-83,  is  given  by  the  gradient  in  Equation  2-12 

E*  (I<3))  =  -(x-i)(4ir)-lv'l(3). 
s  2n  2n 


(2-81) 

i-  e )  is 
(2-82) 

(2-83) 

(2-84) 

(2-85) 

(2-86) 

(2-87) 

(2-88) 
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The  partial  derivative  on  pm  gives  a  radial  component, 

ES  (I<3))  =  -(2/tt)  (*-l)8pS4>  Sz  costa  eikss~2  (ik  -s_1)b  ,  (2-89) 

s  p  2n 

and  differentiation  on  4>m  gives 

E| '  (1^)  =  (2/tt)  (k-1)  8p84>8z  cos^  eikss~2(ik-s*1)pn  sin(4>n-<f>m).  (2-90) 


The  x  component  is 

Es  (I<3))  =  Es  (l(3))cos^n  -  E<3)  (Ph  sin<f>n 
sx  2n  s p  2n  s<f>  2n 


(2-91) 


The  electric  field  from  -  I^2)  is  given  by  applying  the  gradient  operator  to  the 
expression  in  Equation  2-87V  The  first  term,  that  is  the  factor  multiplying  1,  is  identical 


to  E3  (il3))  in  Equation  2-91.  The  second  term  gives  a  radial  component 
sx  2n 

Es  (Ix  -  I2)'  =  -F  eiks  s-3  b  [2bs_1  +  Cnm  +  ikb] 
sp 


(2-92) 


when  Cnm  is  cos (<f>n  -  <^»m),  and  F  is  -2(i/vr)  (x-l)  kpn8z5<f>8p cos<t>n.  The  prime  denotes 
that  the  contribution  of  the  first  term  is  absent.  The  second  term  gives  a  <p  component 

Es4,(i1  "  r2>'  =  'F  eiks  s'2  sin(</>n  -  <t>m)  [l  +  bs“lpn  (ik-s_1)-bpns'2]  (2_93) 

If  we  combine- the  results  we  obtain  for  the  x-component: 

ES  =  v-  i(l) i(2)  - 1(3) 
s  v  2n  2n  2n 


=  [-2  s-1(ik  -2_1)bcos«f>m  -2  s_1  (ik-s-1  )pnsin(<f>n-$m  )sinc  <t> 


m 


+i  s“2kpnb  2bs-1+  cos(ff>n  -  (ffa)  -  ikb]  cos«£m  (2-94) 

+1  s-1kpn  sin(4»n  -  <J>m)  [2b  pn.s_2-l  -ibks-1  Pn]  sin^m]  G* 
where  G  is  [(2 hr)  |k  — X|  8p8<t>  8z  cos$n  e*kss-1] 

Equation  2-94  reduces  to  the  formula  developed  for  the  hollow  cylinder;  where  all  the  pn 
have  equal  magnitude. 


■f,  -  ^  • 
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2.5  EXTERNAL  FIELDS 

The  field  outside  the  dielectric  region  is  found  by  evaluating  Equations  2-26  and  2-40. 
The  integration  extends  over  the  dielectric  region.  It  can  be  done  after  the  total  field 
inside  the  dielectric  is  determined  by  solving  the  integral  equation. 

The  field  near  a  cylinder  is  evaluated  to  test  the  approximations  in  the  derivations.  The 
near  field  also  is  useful  in  understanding  the  effects  of  a  dielectric  cylinder  or  a  radome 
on  an  antenna. 

The  field  from  the  vector  potential  is 

ES  (A)  =  iuA. 


A  measurement  requires  a  probe.  We  consider  a  half-wave  dipole  and  weight  the 
contribution  of  each  cell  by  the  far-field  pattern.  To  simplify  the  expressions,  consider 
the  field  on  the  z-axis.  The  vector  potential  gives  for  the  x  component  (from  one  ring) 

ES  (A^)  *  (k-1)  k“  (4? r)  1  JgE^  cos  cos  9)  (sin  9)  1  pdpd  <p dz  (2-95) 

where  the  factor  cos  (^72  cos  0  )/sin  Q  is  the  far-field  pattern  of  the  probe.  Consider  the 
distance  large  enough  so  that 

r  =  R  +  (acr  -  ap  cos  <p  -  zz')  R-1  (2-96) 


where  a  is  the  mean  radius  of  the  cylinder  and  R  is  (a2  +  z'2)1//2.  Furthermore 
approximate  sin  0 ; 

2  2 

sin  9  =  1-  (a/R)  cos  <p  (2-97) 


2.5.1  FROM  AVERAGE  INTERNAL  FIELD 

To  develop  a  simple  formula  we  omit  the  dependence  of  E  on  $  and  use  an  average  value 
E0.  The  integrations  are  elementary,  so  that 


ES  (A  )  =  2  (k-1)  k"  a  <5p  5z  E  R  *  sine  (ka  6p/R)  sine  (kz1  <5z/R)  P 


■fc’4 


1  -  (a/2R) 

P  is  the  probe  correction. 

The  field  from  the  scalar  potential  is 


(2-98) 


E^  ($)  =  (4*)  * 


(k_1)  ex  dV 


(2-99) 


2-18 


(2-100) 


The  procedures  are  like  those  in  Paragraphs  2.3.2  and  2.4.2.  So 

2 


where 


2  ikR 

E>  a  Eq  f  e1  [. ..]  cos"  <t>  da  dp  dz  , 

R 


(...)  =  (ac  -  ap'  cos  <p  -  zz')  R 


and 


-1. 


[. . .  ]  =  [6  (<r  +  6p)  -  6  (a  -  6p)  -  a  ]  . 

After  doing  the  integrals,  we  have 

ESx  (*x)  =  i  7r  (k-1)  a2  6z  <5p  eikR  (R2\)  ^  (^j  sinc 


(2-101) 


(2-102) 


P'  E 


o 

(2-103) 


Note  the  probe  correction  P'  is  [  1  -  3/8  (tt2/4  -  1)  a2/R2];  it  differs  from  that  in 
Equation  2-98. 


The  total  field  is 


T  I  s  _s 
E  =  E  +  E  (A)  +  E  ($) 

X  X  X 


(2-104) 


The  observable  quantities  are  intensity  |  ET  |  2  and  phase,  the  argument  of  ET. 

X  X 


This  theory  does  not  include  the  effects  of  fields  reflected  from  the  probe  to  the 
scatterer  and  scattered  back  to  the  probe. 

The  formulas  generalize  readily  to  several  rings  forming  a  cylinder.  For  example  for  two 
rings 

TTS  /A  \  -  1 TS  /A  x  /l 7  aikZ'  5Z/R  T7  ikzf  5z  1R\ 

W  -  E  <Al\Eoie  "E02  6  )  (2-105) 


where  Egi  and  Eq2  are  the  fields  in  each  of  the  rings  and  Es  (A^)  is  from  Equation  2-98. 
The  scalar  potential  gives 


Ex  <*2>  ■  E‘  <«2>  (E01 

where  Es  (*2)  is  from  Equation  2-102. 


ikz'  6 z/R  „  3ikz*  5z/R\ 


+  E02e 


/ 


(2-106) 
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2.5.2  FROM  CIRCUMFERENTIALLY  VARYING  INTERNAL  FIELD 

The  external  scattered  field  is  computed  from  a  general  formula,  which  gives 

E^AI  «s*n(Pn3Kn/S*n.gn9Pn/9xn), 

where  A  is  -4  (  *  -  1)  S  z  8$8p/ir,  the  probe  correction  factor  is 

P  °  '  <xn  "  *r)2  rn2]  1/2  [1,<3tn  '  !‘r)/2rn]  ’ 

and  gn  is  g  with  rn  equal  to  [(xn  -  x1)2  +  (yn  -  y1)2  +  (zn  -  z1)2]  1/2. 


(2-107) 
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3.  RESULTS  FROM  SCALAR  GREEN'S  FUNCTION  METHOD 


3.1  CYLINDERS 

3.1.1  INTERNAL  FIELDS 

The  first  example  is  a  single  ring,  diameter  2a  equal  X,  length  2Sz  =  0.16 X,  dielectric 
constant  2.6.  Because  incidence  was  axial,  necessary  conditions  on  the  computed  results 
were  that  the  field  at  0  =  0  equal  that  at  4>  =  tt  and  that  the  field  at  <t>  =  tt/2  equal 
that  at  <i>  =  3ir/2.  This  condition  was  well  satisfied. 

Figure  3-1  shows  computed  internal  values  of  for  several  values  of  N;  the  graph 
suggests  subdivision  into  24  cells  is  adequate  for  stable  values. 

Figure  3-2  shows  computed  internal  values  of  ET  for  a  cylinder,  also  of  diameter  X,  but 
with  length  0.32  X.  Two  rings  were  utilized. 

3.1.2  EXTERNAL  FIELD  FROM  CIRCUMFERENTIALLY  UNIFORM  INTERNAL 
FIELDS 

The  approximation  of  Section  2.5  was  utilized  to  compute  external  fields.  Figure  3-3 
shows  computed  values  of  ET  outside  the  ring,  on  the  z-axis,  with  and  without  probe 


Figure  3-1.  Phase  and  Amplitude  of  ET  for  Single  Ring:  8z  =  0.100  in.,  8p  =  0.030  in., 
x  =  2.6,  frequency  =  9.375  GHz.  For  <f>  =  0°  (o),  <t>  =  90*(*) 
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Figure  3-2.  Phase  and  Amplitude  of  for  Two  Rings.  8z  =  0.1045  in.; 

8p  =  0.0313  in.  at  center  of  second  cell.  Symbols  as  in  Figure  3-2. 

corrections.  Figure  3-4  shows  computed  valued  of  ET  outside  the  cylinder.  Figure  3-5 
shows  computed  values  for  the  cylinder,  but  the  area  of  the  cylinder  cross  section  was 
reduced  by  80%  in  evaluating  the  external  field  integrals,  Equations  2-105  and  2-106. 
The  area  reduction  makes  the  volume  of  the  cylinder  equal  to  that  of  the  elementary 
integration  cells. 
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Figure  3-3.  Total  Field  Behind  Dielectric  Ring.  Length:  0.2  in.;  thickness  26 p: 
0.06  in.;  mean  diameter  2a:  1.25  in.;  dielectric  constant  2.6.  Measured  with 
dipole  probe:  ( • );  computed  without  probe  corrections  (o);  computed  with 
probe  correction  (+).  Internal  field:  1.19  exp  i  7.8*  for  N  =  64 
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Figure  3-5.  As  in  Figure  3-4  but  with  8p  Reduced  by  20  Percent 
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3.1.3  EXTERNAL  FIELDS  FROM  CIRCUMFERENTIALLY  VARYING  INTERNAL 
FIELDS 

This  section  gives  external  field  values  computed  from  Equation  2-107,  which  includes 
the  circumferential  variation  of  the  internal  field.  Figures  3-6  through  3-15  show  the 
intensity  jET|  2  and  phase  A^>  of  the  total  field  on  transverse  paths  at  distances  0.5, 
0.75,  and  1.0  wavelengths  from  the  center  of  a  ring  of  length  28z  =  0.21".  Incidence 
was  axial. 

Figure  3-16  shows  axial  dependence. 

3.1.4  POLARIZATION  DEPENDENCE  OF  INTERNAL  FIELDS 

Figure  3-17  shows  computed  values  of  I ET  |  2  inside  one  or  two  rings.  The  graph  shows 
data  for  a  cell  in  the  E-plane  and  for  a  cell  in  the  H-plane.  The  E-plane  cells  are  at  the 
sides  of  the  cylinder,  and  the  H-plane  cells  are  at  the  top  and  bottom;  recall  the  field 
was  horizontally  polarized.  Clearly  a  polarization  dependence  exists. 
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Figure  3-6.  Total  Field  Intensity  for  a  Ring;  z  =  0.5 X;  2 5p  =  0.625", 
25z  =  0.21",  p  =  0.594";  *  =  2.6.  (-):  Computed;  (X):  Measured 
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Figure  3-7.  Total  Field  Phase  for  a  Ring  like  that  in  Figure  3-6. 
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Figure  3-11.  As  in  Figure  3-7  but  z  =  A. 
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Figure  3-12.  Total  Field  Phase  for  a  Cylinder  Composed  of  Two  Rings; 
Length  0.42".  (-):  Computed;  (X):  Measured. 
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Figure  3-13.  Total  Field  Intensity  for  a  Cylinder  Composed  of  Two  Rir 
=  0.21";  28p  =  0.0625;  P  =  0.594";  «  =  2.6.  (-):  Computed;  (X):  M 
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Figure  3-14.  As  in  3-12,  but  z  =  X . 
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Figure  3-15.  As  in  Figure  3-13,  but  z  =  A. 
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Figure  3-16.  Total  Field  Behind  Dielectric  Ring  with  Dimensions  Given  in  Figure  3-6 
Measured  (• );  Computed  with  Equation  2-107(+). 
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Figure  3-17.  Computed  Internal  Values  of  |Et|  2  for  1  and  2  Rings.  28p  =  0.0625", 
2Sz  =  0.21",  p  =  0.594";  N  =  64  Per  Ring. 

3.2  CONE 

Calculations  were  done  for  the  cone  in  Figure  2-4;  dimensions  are  given  in  Figure  2-5. 
Incidence  was  axial;  frequency  was  9.400  GHz;  polarization  was  horizontal  (parallel  to 
the  x-axis  of  Figure  2-3,  with  the  cone  axis  parallel  to  the  z-axis.  The  cone  had  13  rings 
although  the  cone  in  Figure  2-4  has  15. 

Calculations  were  first  done  for  11  rings  as  well  as  13,  with  the  subdivision  given  in 
Table  3-1,  and  for  thickness  of  all  rings  2Sp  =  0.63". 

Computed  values  of  total  field  intensity  within  the  cone  are  shown  in  Figure  3-18. 

The  computed  data  in  Figure  3-18  are  for  two  generators.  One  is  the  E-plane;  the  other 
the  H-plane.  The  polarization  dependence  is  large,  especially  near  the  tip.  The  nature 


Table  3-1.  Computational  Parameters 


NUMBER  OF  RINGS 

CELLS  PER  RING 

13 

12 

13 

8  in  smaller  6  rings,  16  in  larger  7  rings 

11 

16 

CDNICRL  SHELL  1 NTERNRL  FIELD 


Figure  3-18.  Computed  |E^|  2  inside  conical  shell.  The  inset  shows  the  loci 
of  points.  E  refers  to  points  in  the  x  -  z  plane,  the  E-plane;  H  refers  to 

y  -  z  plane,  the  H-plane. 
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of  the  dependence  resembles  that  in  Figure  3-17  for  two  rings  of  equal  radii.  The  curves 
in  Figure  3-18  show  that  the  numerical  representation  influences  |  ET  |  2.  However, 
these  differences  are  smaller  than  those  between  the  curves  for  the  E  and  H  planes. 
Near  the  tip,  in  the  region  z  <  -1,  the  intensity  is  -10  dB  below  the  incident  intensity, 
but  the  H-plane  intensity  is  4  dB  to  5  dB  above  the  incident.  Near  the  base  the  intensity 
exceeds  the  incident  for  both  planes. 

The  H-plane  curve  for  8  cells  in  the  smaller  rings  differs  from  the  curves  for  the  larger 
number  of  cells  per  ring. 

The  fact  that  |  ET  |  2  is  greater  in  the  H-plane  than  the  E-plane  is  reasonable  because 
the  electric  field  is  approximately  parallel  to  a  cell  in  the  H-plane  but  approximately 
orthogonal  to  a  cell  in  the  E-plane.  However  the  large  dependence  on  polarization  is 
smaller  for  a  single  ring. 

The  smallest  rings  have  thickness  approximately  that  of  the  radii,  violating  an  assump¬ 
tion.  Therefore  calculations  were  made  by  representing  each  of  the  three  smallest  rings 
by  a  pair  of  concentric  rings.  Figure  3-19  shows  the  computed  internal  field  intensity 
|  Et|*.  The  values  for  the  E-plane  depend  on  radial  position,  but  those  for  H-plane 
depend  negligibly  on  radial  position.  The  results  in  Figure  3-19  also  differ  from  those  in 
Figure  3-18  near  the  base. 

The  external  field  was  computed  for  a  transverse  path  at  z  =  1.812",  which  is  a  half 
wavelength  behind  the  base  of  the  cone.  Figure  3-20  shows  the  results.  The  subdivision 
of  the  three  rings  nearest  the  tip  influences  the  field. 
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Figure  3-20.  | ET| 2  at  Z  One-Half  Wavelength  Behind  Largest  Ring 
[z  =  1.812").  Measured  (X);  Computed  ( •)  for  13  Rings,  12  Cells/Ring 
5p  =  .065";  Computed  with  16  Rings,  Smallest  Rings  Have  8p  =  .033"  (•). 
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4.  MEASUREMENT 


Intensity  and  phase  near  the  cylinders  and  cone  were  measured  with  a  setup  that  included 
a  half-wave  dipole  and  a  network  analyzer.  Measurements  were  made  for  horizontal 
polarization  in  the  system  of  Figure  2-1.  The  probe  was  scanned  in  the  x  direction  or  y 
direction  for  fixed  values  of  z. 

For  the  cylinders  Figures  3-3  through  3-6  show  two  observed  values  on  the  z  axis;  one 
value  comes  from  each  scan.  Figures  3-6  through  3-15  show  the  data  on  transverse 
scans. 

For  the  cone,  Figure  3-20  shows  measured  values  of  phase  and  intensity  for  z  =  1.816" 
as  a  function  of  x.  This  scan  is  half  a  wavelength  behind  that  base  of  the  cone. 
Figure  4-1  shows  intensity  on  two  orthogonal  scans  at  the  center  of  the  second  largest 
ring. 
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Figure  4-1.  Phase  and  Intensity  Measured  in  Plane  of  12th  Ring,  z  =  0.86"  for 
13-ring  cone.  The  electric  field  was  orthogonal  to  the  y-z  plane  of  Figure  1. 
The  H  plane  is  the  y-z  plane  (x  =  0)  and  the  E  plane  is  orthogonal  to  the 
H  plane  (that  is,  y  =  0).  The  probe  was  a  half-wave  dipole. 


5.  DISCUSSION 


5.1  CYLINDERS 

Figure  3-1  suggests  that  convergence  is  obtained  for  a  single  ring  divided  into  twenty- 
four  cells,  but  Figure  3-2  suggests  that  two  contiguous  rings  may  require  a  somewhat 
larger  number  of  cells  per  ring.  The  phase  in  the  H-plane  generator  is  the  most  sensitive 
quantity. 

The  axial  dependence  in  Figure  3-3  shows  the  phase  of  the  external  field  is  somewhat 
overestimated.  Discrepancies  are  28  percent  at  the  closest  spacing  of  probe  and  ring  but 
are  a  few  percent  at  wavelength  spacing.  Probe  correction,  in  the  theory,  reduces 
discrepancies.  Figure  3-5  shows  that  a  reduction  of  Sp  by  20  percent  reduced  discrep¬ 
ancies  to  small  values,  comparable  to  the  measurement  precision.  The  reduction  of  Sp  is 
possibly  related  to  the  approximation  made  in  evaluating  the  integral  in  Paragraph  2.3.1, 
specifically  to  the  integral  over  8 . 

The  axial  dependence  in  Figure  3-16  is  comparable  to  that  in  Figure  3-3  despite  the 
approximation  of  using  circumferentially  independent  fields  for  the  results  shown  in 
Figure  3-3. 

The  transverse  dependence  in  Figure  3-6  shows  discrepancies  of  five  percent  in  |ET|  at 
off-axis  positions  of  approximately  one  inch.  In  general,  phase  discrepancies  are 
approximately  20  percent  for  small  probe-ring  separation.  The  shape  of  the  measured 
and  computed  curves  agree  well. 

5.2  CONE 

Figure  3-18  shows  a  polarization  dependence  of  the  internal  field  I  E^l .  The  approxima¬ 
tion  that  required  ring  radius  to  be  much  larger  than  thickness  was  tested.  Figures  3-18 
and  3-19  show  that  this  approximation  produces  significant  internal  field  changes. 

The  measured  data  in  Figure  4-1  show  a  polarization  dependence.  The  H-plane  intensity 
exceeds  the  E-plane  intensity  by  approximately  1-1/2  dB.  This  difference  is  approxi¬ 
mately  that  for  the  internal  fields  of  the  larger  rings;  see  Figure  3-19.  Figure  3-20 
shows  that  discrepancies  are  reduced  by  performing  calculations  with  smaller  ring 
thicknesses,  for  the  rings  near  the  tip  of  the  cone.  Tlie  differences  between  measured 
and  computed  phase  are  large  for  the  probe  on  the  axis,  but  fall  to  smaller  values  at  the 
first  minimum.  The  shape  of  the  measured  and  computed  curves  agree,  especially  for 
the  graph  computed  with  small  values  of  Sp  near  the  tip.  It  is  unclear  whether  the 
discrepancies  are  caused  by  numercial  or  theoretical  approximations. 
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6.  TENSOR  GREEN'S  FUNCTION  FORMULATION 


6.1  INTEGRAL  EQUATION 

Maxwell's  equations  for  time  harmonic  electric  and  magnetic  fields  of  angular  frequency 
a)  in  the  presence  of  a  nonmagnetic  scatterer  of  dielectric  function  k  (r)  can  be  combined 
to  yield  a  wave  equation  for  the  electric  field  E(r): 

VxVx£  -  k2E  =  k2  Oc-l)£  (6_1) 

Here  k  =  2ir/\,  =oVc  is  the  free  space  wavenumber,  X0  the  free  space  wavelength,  and  c 
the  speed  of  light  in  vacuum.  The  tensor  (or  dyadic)  Green's  function11!’: 

F(k,  r' ,  r)  =  (£  +  k-2  V  V)  (exp  (ik  |  r'  -_r|)  (47r|r'  -_r|f 1  (6_2) 

satisfying 

VxVx  T-k2  r  =  1  /(r-rl 

V.£  =  0  (£  *  _r') 


3 

(where  6  (r  -  r')  is  the  three-dimensional  Dirac  delta  function)  can  be  used  to 
immediately  rewrite  (6-1)  as  an  integral  equation: 


E  (k,  r')  -  E1  (k,  r')  +  k2  /  dV  T(k,  r',  r)  E  (k,  r)  [k  (r)  -  l]  (6-3) 

V 

s 

Here  E*  is  the  incident  electric  field  satisfying  the  homogeneous  free  space  equation 


I  2  I 

VxVxE  -  k  E  =  0 

and  the  integration  extends  over  the  volume  V  of  the  scatterer  (everywhere  x(r)  *  1). 

To  begin  the  transition  to  the  moment  method  approximation  of  (6-3),  we  partition  V 
into  i?  subvolumes  ("cells")  each  of  volume  V^;  moreover,  we  assume  «(r)  =  <c  =  constant 
inside  the  scatterer.  Then  (6-3)  becomes: 


£(k,  r)  =  E1  (k,  r')  +  k2 


fc-1)  J  dV  X  (k,  r',  r)  •  E  Oc,  r) 


n  * 


n  ,V 
r  n 


+k2  (k-1)  f  dVr  (k,r',r)  £(k,r) 


nr' 


(6-4) 


6-1 


where  Vn  ,  is  the  volume  of  the  cell  containing  the  observation  pointer".  If  H  is  outside 
the  scatterer,  the  last  integral  in  (6-4)  is  absent. 

Now  we  introduce  the  essential  assumption  underlying  the  moment  method:  the  nth  cell 
is  assumed  to  be  small  enough  that  the  electric  field  E(r)  may  be  approximated  as  a 
constant  jgn  throughout  the  cell.  As  a  consequence  E(rO  may  be  removed  from  the 
integral  in  (6-4): 


g  (k,  r')  =  E1  (k,  rO  +  k2  (*-l)  E  •  /dV  £  (k,  r',  r) 

/  -y  n  V  ~~ 


n  *  n  . 

r 


(6-5) 


+  k2$c-l)£  •  / dV  (k,r',r) 

n  ,, 


nr' 


The  major  remaining  task  is  the  evaluation  of  the  two  integrals  appearing  in  (6-5): 


£(k,  r',  r)sf dV  £  (k,  r',  r)  (6-6) 

V 

n 


Substituting  (6-2)  for  £  into  (6-6): 


I  (k,  r',  r)  =  (4fff 1  / dV  (1  +  k~2  v'  VO  g  (k,  r',  r)  (6-7) 

*  vn 

where  g(k,  r'  r)  =  (|rf  -£|)  1  exp  (ik  )r_'  -_rj)  is  the  scalar  Green's  function.  To  actually 
perform  the  integration  (6-7)  we  approximate  the  nth  cell  as  a  sphere  of  radius  a 
centered  at  r  .  If  we  temporarily  exclude  the  case  n  =  n^.  so  that  the  observation  point 
£  is  not  within  the  cell,  we  may  remove  the  gradient  operators  from  the  integral: 


Hk,  r',  n)  =  (4*rf 1  (i  +  k-2  vV)  /  dVg(k,  r',  r) 


The  integration  over  the  spherical  volume  is  conveniently  carried  out  by  Fourier 
decomposing  g;  the  result  is 


£(k,  r',  =  k~3  (sink  a  -  kacoska)  (£  +  k-2  v'v')  g  (k,  r'.  r  )  (lx'  -x  1 2a)  (6-8) 


6-2 


When  n  =  nr>,  the  source  point  £  and  observation  point  r'  coincide  in  the  integration 
volume  Vn^and  the  integrand  in  (6-7)  is  formally  divergent.  However,  a  principal  value 
prescription1^  may  be  applied  to  give  (6-7)  a  unique  value  in  this  case: 


Hk,  r',  nr,)  =  -k  2  (ika  -  1)  exp  (ika)  +  lj  1_ 


In  view  of  (6-6),  (6-8),  and  (6-9),  equation  (6-5)  becomes: 


E  (k,  r')  '  E*  (k,  r)  +  1)  k  1  (sinka  -  kacosca) 


ZE  •  (1  +  k'2  v'vO  g  (k,  r',  r  ) 
— n  ~  n 


n  *  n  , 
r 


(6-9) 


(6-10) 


+  (1-x)  f|  (ika  -  1)  exp  (ika)  +  lj  En 

The  last  term  is  absent  if  r'  is  outside  the  scatterer.  Rearrangement  of  (6-10)  yields 
E  (k,  r')  a  [|  (x-  l)  (ika  -  1)  exp  (ika)  +  xj-1  E1  (k,  r') 


+  (x-  1)  k-1  (sinka  -  kacoska)  (x-  1)  (ika  -  1)  exp  (ika)  +  xj 


(6-11) 


:E  V  %  + 

n*  n 


k~2  vV)  g  (k,  r',  rR) 


(valid  for  r'  within  the  scatterer) 


E  (k,  r')  =  E1  (k,  r')  +  (x-  1)  k  1  (sinka  -  kacoska) 


x£  En  •  a  +  k-2  vV)  g  (k.  r'.  rn> 


(valid  for  £"  external  to  the  scatterer) 


(6-12) 
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From  (6-11)  we  obtain  for  the  field  in  the  mth  cell  in  the  scatter: 


~m  ~  [I  ^ika  ~  ^  exp  ^ika^+  K]  1  |-^i  *  k  1  (sinka  ”  kacoska) 


B*- 


1)  ( ika  -  1)  exp  (ika)  +  * 


]“«Ev 


(1  +  k  Vm  v  )  g  (k,  R  ) 
m  m  mn 


(6-13) 


n  *  m 


where 


*k-Rnn)  *  Rmn  exp  <ik  Rmn> 


Rmn  '~mJ  ~m  ~n' 


Since 


k'2Vm  7m  «  (k’  R  „)  =  g  (k,  R__)  jk'2  R“fn  [3  (1  -  ikR  )  -k2P2  1  R  £r 
v  m  ▼  m 0  rnn  mn  |  mn  L  mn  mnj  mn  r 


+  (kR  )~2  (ikR  -1)  1  (R  *  Q) 
mn  mn  ^  |  mn 


Equation  (6-13)  becomes 


~m  ~  £|  (k-i)  (ika  -  1)  exp  (ika)  +  1  E ^ 

+  («-l)  k  3  (sinka  -  kacoska)  [|  (K“l)  (ika  -  1)  exp  (ika)  +  k]-1 

x^  g(k,  R  )  R-2  |k2R2  +  ikR  -  ll  E 

-f;i  mn  (  mn  l  mn  mn  J  — n 


x  >  .  8 <k> R  > 


12  2 

kR  +  ikR 
1  mn  mn 


-iK 


(6-14) 


-4  r  22 

+  R  3  (1  -  ikR  )  -k  R 
mn  L  mn  nr 


) -k2R2  1  (R  •  E  )  R 
mn  mnj  ~mn  ~n  —  mn  | 
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Equation  (6-14)  is  our  main  result.  Note  that  for  ka  «1: 


j  fc-1)  (ika  -  1)  exp  (ika)  +  <  ~  (<+2)/3 


-3  3 

k  (sinka  -  kacoska)  =:  a  /3 


so  that  in  this  limit  (6-14)  is: 

=  3  (<c+2)“^i  +  (k+2)'1  (k-1)  a3 


•m 


xY]  g  (k.  R  )  j  R~2  [*k2R2  +  ik  R_1  1  E 

mn  (  mn  L  mn  mnj  ~n 


(6-15) 


R'4  1*3  (1  -  ik  R  )  -  k2R2 
mn  L 


mn 


J (s- 


...  E  )  R  }  (ka«l) 
mnj  ~mn  -n  "mn  * 


Now  consider  a  collection  of  oscillating  electric  dipoles.  The  electric  field  at  £  due  to 
one  such  dipole  is11: 

^Dipole  _  k2  |  x  pj  x  (r/p) |  g  (k,  r) 

+  1 3  (r/r)  [(r/r)  •  P]  -  P  j  (r  2  -  ikr  X)  g  (k,  r) 

=  g  (k,  r)  |r  2  |j<2r2  +  ikr  -  l]  P  +  r  4  [3  (1  -  ikr)  -k2r2J  (r  •  P)_r  j 

where  P  is  the  dipole  moment.  If  the  dipole  moment  is  induced  by  an  electric  field 
jrinternal^  we  have 

p  _  .^internal 

where  Y  is  the  "molecular"  polarizability  (i.e.,  the  polarizability  of  the  dipole  when 
immersed  in  the  collection  of  all  the  other  dipoles).  By  the  Claussius-Mossotti 
equation14: 

Y  =  (4  ir  N)"1  3  (K-l)  (K+2)"1 

where  N  is  the  number  of  dipoles  per  unit  volume.  Assuming  each  dipole  occupies  a 
sphere  of  radius  a, 


implying 


N  =  (4  a3/3)-1 


Y  =  (k-2)'1  (k+1)  a3 
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Therefore, 


^Dipole  .  (|(+2)-l  a3  g  {kj  r)  j  r-2  j^2p2  *  ikr  .  ,]  g> 

i>r  | 


internal 


(6-16) 


1  I  internal 

Further,  if  E  is  the  incident  field  in  free  space,  then  the  internal  field  G  in  any 


dielectric  sphere  due  to  E*  is1®: 


EI  internal  =  3  ((C+2)-l  £I 


(6-17) 


Comparing  (6-16)  and  (6-17)  to  (6-15),  we  observe  that  for  ka«l,  the  moment 
approximation  consists  of  replacing  a  dielectric  scatterer^by  a  collection  of  interacting 
electric  dipoles  each  occupying  a  spherical  volume  Air  a  /3  and  whose  induced  dipole 
moments  obey  the  Clausius-Mossotti  equation.  This  result  is  the  same  as  that  of  Purcell 
and  Pennypacker16  (with  the  exception  of  the  incident  field  term,  where  they  assume 
gl  inlernal  =  E1),  and  so  constitutes  an  independent  verification  of  our  main  result  (6-14) 
in  the  long  wavelength  limit. 

6.2  COMPUTATIONAL  TECHNIQUE 


Eq.  6-14  may  be  written  more  compactly  as 


E  =  6_1EI  +  <5“10! 
~m  -m 


where 


n  #  m 


_3 

a  =  (k-1)  k  (sinka  -  kacoska) 


V  jfl  E  +  Ymn  •  E 
/  ,  (pmrr~n  =r  ~-n 


(6-18) 


smn  *  S  <"•  Rmn>  Rmn  [k2RL  MkRmn -1] 

Rm„>  Rmn  [3  11 '  “J  5, 


(k-i)  (ika  -  1)  exp  (ika)  +  « 

Note  that  /3  and  7  are  symmetric  ij  by  17  matrices. 
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Introduce  i,  j  for  3-vector  indices.  In  component  form  (6-18)  becomes: 


E  .  =  +  * _1a  „  -J3_ _  *  +  vmnii5  Enj 


'  £  £^mn*ij  +  Vmnii  ni 
n  *m  j 


By  defining 


Ln  ,  V(1“6mn) 


9  ..  =  y  ..  (l  -  6  ) 

mmj  mni]  mn 


we  may  formally  extend  the  sum  in  (6-19)  to  include  n  =  m: 


E  .  =  <5"1EI  .  +  <5-1  or  (£  9  ..)  E  . 

mi  mi  2j  L  mmj  mmj  nj 

n  j 


r?  3 

Emi  ~  H  [6firnn  5ij  "  ^mn  6  ij  +  Ymnij^]  En 

j  » 1 

If  we  introduce  "composite"  indices  M,  N  by  the  mapping 


M  =  3  (m  -  1)  +  i 


N  =  3  (n  -  1)  +  j 


(1  <  M,  N  <  3T?) 


Eq.  (6-20)  takes  the  form 


T  0TJ 

€  =  T  A  € 

M  v/-  MN  N 
N  =  1 


where 


'M  =  6  3  (m  -  1)  +  i 


(6-19) 


(6-20) 


(6-21) 


il  a  E1 . 
J  tm 


SN  *  ,[3<n-1)*i]‘Enj 

AMN  -  4  (m  -  »  *  i]  b  <n  -  h]  -  S5MN  6y  -  V  6  y  ‘  \mij) 


Amn  is  a  symmetric  3tj  by  3t?  matrix. 

According  to  (6-21),  the  moment  method  approximation  has  reduced  the  problem  of 
finding  the  electric  field  in  a  dielectric  scatterer  to  the  algebraic  task  of  solving  3 17 
simultaneous  equations  for  317  complex  unknowns.  A  variety  of  practical  computational 
techniques  exist  for  solving  such  a  system  of  equations  for  17  as  large  as  103  -  104  cells. 

The  only  remaining  latitude  in  the  approximation  is  the  selection  of  a,  the  radius  of  the 
spherical  cells,  and  1?,  the  total  number  of  such  cells.  Since  we  require  that  £  be  almost 
constant  in  any  one  cell,  we  make  this  requirement  quantitative  by  asking  that  the  phase 
change  in  E  across  any  cell  diameter  2a  be  less  than  .2  radians.  Thus: 

2a /U/Jk)  <.2 


or 


a  < 


(6-22) 


If  we  further  assume  that  the  total  volume  of  the  tj  spherical  cells  equals  the  volume  V 
of  the  scatterer: 


V  =  V  /&  rra3)  >  V  (10^)3/(|  n\  3) 
S3  —  s  3  0 


(6-23) 


7.  RESULTS  FROM  THE  TENSOR  GREEN'S  FUNCTION  FORMULATION 

An  earlier  report1**  treated  scattering  from  a  dielectric  rod  and  dielectric  ring  utilizing 
the  scalar  Green's  function  technique.  In  this  section,  we  reexamine  those  cases  using 
the  Green's  tensor  formulation  developed  in  Section  6  of  this  report.  The  measurements 
shown  below  are  taken  from  reference  10,  where  they  are  discussed  in  more  detaiL 

7.1  ROD 

The  first  example  considered  was  a  dielectric  rod  (/c  =  2.6)  of  length  1.9  inches  and 
diameter  0.139  inches  aligned  lengthwise  along  the  z-axis.  Figure  7-1  shows  the  rod  and 
its  approximation  as  15  overlapping  spheres,  each  of  radius  .0765  inch.  The  incident 
wavelength  was  1.259  inches  and  both  parallel  and  perpendicular  incident  polarizations 
were  considered.  Figures  7-2  and  7-3  display  the  computed  internal  fields  for  the  two 
polarizations.  For  comparison,  note  that  the  internal  field  in  an  infinite  rod  of  the  same 
diameter  and  dielectric  constant  is  1.4  dB  above  the  incident  field  in  the  case  of  parallel 
polarization.  Figures  7-4  through  7-11  exhibit  the  calculated  and  measured  external 
fields  at  two  different  distances  from  the  rod.  With  the  exception  of  Figure  7-9,  the 
agreement  between  theory  and  measurement  is  good.  Since  we  would  expect  the 
measured  phase  difference  to  vanish  as  z  becomes  large,  the  discrepancy  in  Figure  7-9  is 
evidently  a  systematic  measurement  error. 
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Figure  7-1.  Dielectric  Rod  as  Approximated  by  15  Dielectric  Spheres. 
Rod  Length:  1.9",  Rod  Diameter:  0.139";  Sphere  Radii:  0.0765". 
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Figure  7-2.  Computed  Internal  Field  Intensity  for  the  Dielectric  Rod  of  Figure  7-1. 
(O):  E1  Parallel  to  Rod;  (*):  E1  Perpendicular  to  Rod. 
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Figure  7-3.  Computed  Internal  Field  Phase  for  the  Dielectric  Rod  of  Figure  7-1. 
(O):  E*  Parallel  to  Rod;  (*)  E*  Perpendicular  to  Rod 
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Figure  7-4.  Total  Field  Intensity  for  the  Dielectric  Rod  of  Figure  7-1  for  E 
Parallel  to  Rod  and  x  =  .5A.  (-):  Computed;  (X):  Measured. 
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Figure  7-5.  Total  Field  Phase  for  the  Dielectric  Rod  of  Figure  7-1  for  E 
Parallel  to  Rod  and  x  =  .5A.  (-):  Computed;  (X):  Measured. 
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Figure  7-6.  As  in  Figure  7-4,  but  x  =  1.5X. 
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Figure  7-7.  As  in  Figure  7-5,  but  x  =  1.5X. 
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7.2  RING 

The  second  example  was  a  dielectric  ring  (k  =  2.6)  lying  in  the  x-y  plane,  as  shown  in 
Figure  2-1.  The  incident  wavelength  was  again  1.259".  Figure  7-12  shows  the  ring 
approximated  as  36  overlapping  spheres,  each  having  radius  0.0684".  Figures  7-13  and 
7-14  depict  the  computed  internal  fields  E||  and  Ex  as  functions  of  angle  around  the  ring. 
Figures  7-15  and  7-16  display  the  computed  and  measured  external  field  along  the  ring 
axis.  The  agreement  is  quite  good  for  z>A  ,  with  small  discrepancies  for  z<  A.  However, 
it  should  be  pointed  out  that  no  attempt  has  been  made  in  the  calculations  to  compensate 
for  the  non-isotropic  receiving  pattern  of  the  field  probe.  Thus,  as  the  probe  moves 
closer  to  the  scatterer,  the  solid  angle  subtended  by  the  scatterer  at  the  probe  grows 
larger,  and  the  effect  due  to  the  non-isotropy  of  the  probe  becomes  more  and  more 
pronounced.  This  non-isotropy  in  combination  with  the  field  disturbances  caused  by  the 
mere  existence  of  the  probe  probably  suffices  to  account  for  the  near  field  discrepancies 
exhibited  in  Figures  7-15  and  7-16. 


Figure  7-12.  Dielectric  Ring  (as  in  Figure  2-1).  Approximated  by  36  Spheres. 

Sphere  radii:  0.0684". 
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Figure  7-13.  Computed  Internal  Field  Intensity  as  a  Function  of  Azimuth  for  the 
Dielectric  Ring  of  Figure  7-12.  (O):  Component  of  ET  Parallel  to  E1; 

(*):  Component  of  ,gT  Perpendicular  to  J21. 


Figure  7-14.  Computed  Internal  Field  Phase  as  a  Function  of  Azimuth  for  the 
Dielectric  Ring  of  Figure  7-12.  (O):  Component  of  ET  Parallel  to  Jg1; 
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APPENDIX  A 

EQUIVALENCE  OF  TENSOR  AND  SCALAR  FORMULATIONS 


In  Equation  2-19  use  Equations  2-9  and  2-10  so  that 


4  ff  ES  s  /  ^iwp^gJ  +  PVfcj  dV  <A— 1) 

The  equation  of  continuity  gives  p  =  (io>)  1  V.Jj  thus 

4ffto«oES  =  i  5  [K^gJ-(  V.  J)vg]  dV,  (A-2) 

2  2 

where  kq  is  <a  fiQ  «Q,  and  observe  the  vg  =  “  Vg- 

The  second  term  in  Equation  A-2  is  17 

(  V  .  J)vg  =Z‘a  V‘(J  -  (J.  V  )  Vg  (A-3) 

a 

In  the  volume  integral  of  Equation  A-l  the  first  term  of  Equation  A-3  can  be  changed  to 
surface  integral.  The  surface  integral  vanishes  outside  a  sphere  that  contains  the 
sources.  Therefore,  Equation  A-2  yields. 

ES  =  i  (4ffft»€0)'1  /  [k^gJ  + J.  VVgjdV  (A-4) 

This  expression  is  given  in  several  references  (see  e.g.,  Reference  9). 
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